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We describe a switching device based on Rabi oscillations. 
The device consists of a well region separated from a free 
region by a potential barrier. The potential landscape is ad- 
justed so that one bound state and one quasi level are present. 
By applying a microwave field with a driving frequency close 
to the separation between the two levels, a particle initially 
in the ground state can be activated to the quasi level, and 
subsequently tunnel to the free region. The probability for 
tunneling in the free region exhibits a plateau structure, as 
the wave function is emitted by bursts after each Rabi oscilla- 
tions. The leakage current can then be controlled by varying 
the amplitude of the external field, the barrier height/width, 
and the frequency mismatch. 



uj r = \I6uj z H ^- 
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The recent advances in the fabrication of semicon- 
ductor nano-structures have generated a rebirth of one- 
dimcnsional quantum mechanics. Indeed, while one di- 
mensional (1 D) scattering problems were considered on 
a purely academical level in the sixties, one now has the 
means to fabricate and design such 1 D structures using 
molecular beam epitaxy techniques, and test the predic- 
tion of relatively simple theoretical models. The motiva- 
tion to study such artificial structures comes to a large 
extent from the hope that these could be ultimately used 
as elementary building blocks in the next generation of 
electronic devices and computers. Unfortunately, this re- 
search is still at its preliminary stage, and the conceptu- 
alization of even the simplest devices already represents 
a considerable challenge. In this letter, we do a step in 
this direction, and we consider a device which allows to 
control the switching of an electric current: the device 
exploits the Rabi oscillations between a bound state and 
quasi levels in the presence of a microwave field. We 
will restrict the discussion to the single electron case: 
electron-electron interactions will be ignored throughout 
the paper. 

Consider a two state system, with ground state en- 
ergy Eg and excited state E\ in the presence of a har- 
monic perturbation. If the frequency of the perturba- 
tion matches roughly the spacing between the two lev- 
els, the system undergoes oscillations with a frequency 
which is much smaller than the excitation frequency u>. 
This Rabi frequency depends on the mismatch 8uo = 
(Ei — Eq)/% — ui between the level spacing and the exci- 
tation frequency, and the amplitude, or matrix element 
Fig of the perturbation ||: 



If we start with the system initially in the ground state, 
transitions to the excited state will occur, but the system 
will ultimately return in the ground state after a period 
T R = 2tt/uj r . 

We propose to use these oscillations in the following 
device. Consider the potential landscape depicted in Fig. 
[I]: A potential well with depth V is separated from a 
continuum region by a thin barrier with height W. We 
adjust the well width a so that there is only one bound 
state in the well. The continuum states with E > 
extend from the well to the free region. If the barrier 
was infinite, we would have a discrete set of states in the 
well. For a finite potential barrier, these states become 
quasi-levels, which have a finite lifetime in the well. We 
choose the barrier parameters in such a way that only 
one quasi level exists for the interval < E < W . By 
applying a dipole field: 



(x, t) = ex cos(cjt) 
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with a frequency which is close to the transition between 
the ground state and the quasi level, the system will un- 
dergo Rabi oscillations (as long as the wave function has 
not leaked totally out of the well region). By adjusting 
the amplitude e of the external field, the excitation fre- 
quency u) and the barrier parameters, we can control how 
much of the wave function penetrates the free region. 

We perform a numerical calculation to probe this de- 
vice. In a first step, we need to determine the bound and 
continuum states for the potential of Fig. [j]. The "free" 
region is chosen to have a width c which is much larger 
than the well and barrier parameters. The wave func- 
tions in all three regions are determined using the con- 
nection formulas for square potential steps, and the corre- 
sponding energies for the bound and "continuum" states 
are found numerically. To identify the energy around 
which the quasi level is centered, we calculate the in- 
tegrated density of all states with E > 0, and select 
the levels for which this quantity is a maximum. Al- 
ternatively, we calculated all matrix elements of the po- 
sition operator between the ground and excited states, 
and select the level for which the probability of transi- 
tion is a maximum (these two procedures give the same 
result). Once the spacing between the quasi-level en- 
ergy Eq and the ground state energy Eg is known, we 
choose the driving frequency uj with the desired mismatch 
Slu = {E Q - E G )/h-u. 
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For the time evolution, the dipole potential of Eq. (||) 
requires special care. Because of the presence of the po- 
sition operator in Eq. (Q), all terms in the Hamiltonian 
do not commute with each other, and the time ordering 
in the evolution operator |3| has to be considered. How- 
ever, by performing a gauge transformation such that the 
scalar potential in the new gauge is zero and the driving 
field is fully included in the vector potential, we are able 
to obtain a Hamiltonian for which all terms commute 
(in each constant region of the static potential V(x)): 
the corresponding vector potential is constant in space, 
and therefore commutes with the momentum operator. 
For an elementary time step, the time evolution operator 
is then written in the Caley form, and a generalization 
of the finite difference scheme of Goldberg, Schey and 
Schwartz Q for the case of a vector potential is derived. 
More details on the method of solution are provided else- 
where (|]. At t = 0, the particle is taken to be in the 
ground state, and at each time step, we compute the in- 
tegrated density in the free region: 

Pi(t) = dx p(x,t) , (3) 

Ja+b 

as well as the overlap | < G\ip(t) > 2 with the ground 
state. 

In Fig. ||a we plot pi (t) for several values of the driv- 
ing field amplitude. The barrier width b = 2 and height 
W = 3 are chosen to be large enough that the "escape 
time" of the wave function is large compared to other 
characteristic times of the problem. Moreover, we have 
chosen the driving frequency to be close to the resonance 
condition {8lo/uo — 0.0001), which is smaller than the 
level spacing for "continuum" energies E > 0, in order to 
check agreement with the two level approximation. For 
the above parameters, the external frequency is roughly 
u> = 1.5 in dimensionless units. The infinitesimal time 
step St for the numerical evolution has to be chosen to 
be small compared to the period of the external field: 
here we choose St = 0.0125 x (2tt/lu). The integrated 
density exhibits steps or plateaux, which allow to iden- 
tify the Rabi frequency lur. Superposed to the plateaux 
structure are small amplitude oscillations which period 
corresponds to the driving frequency lj. As we start the 
simulation, transitions to the quasi-levels and neighbor- 
ing levels start occurring, but after a period Tr = 27r/w^, 
the contribution of the wave function which remained 
trapped in the well has returned for the most part in 
the ground state. This is illustrated in Fig. |2|b : indeed, 
aside from a slow decay associated with the transparency 
of the barrier, | < Ec^it) > | 2 oscillates with period 
which matches exactly the plateaux structure. After the 
first plateau, Pi(t) picks up again as the next oscillation 
returns the trapped wave function to the excited states. 
The wave function in therefore emitted by "bursts" from 
the well region every time a Rabi oscillation has been 



completed. 

We note that upon doubling the coupling strength, the 
period of the oscillations is halved, confirming the fact 
that the Rabi oscillation frequency scales linearly with 
the field amplitude if the resonance condition is met. We 
compare the measured Rabi frequency with the predic- 
tion for the two level case Hujr = \ < G\x\l > \e, on reso- 
nance. To estimate the magnitude of the matrix element, 
we choose the ground and first excited states of an infi- 
nite well, which yields < G\x\l >= 8a/97r 2 . With a = 4, 
the predicted Rabi frequency is therefore uj v r ~ 0.04 (di- 
mensionless units). The measured Rabi frequency yields 
W J? ~ 0.1, which is in reasonable agreement with the two 
level result: the discrepancy between the two arises from 
the matrix element estimate. In addition to these effects, 
we notice in Fig. |]a that after a Rabi oscillation, the in- 
tegrated density is, as expected, larger for higher values 
of the driving field amplitude. 

Next, in Fig. ||a, we compare the evolution of this sys- 
tem for several values of the barrier widths, keeping the 
other parameters fixed. By changing this parameter, we 
effectively vary the escape time of the particle trapped 
in the well. For small widths (b < 1), the escape time of 
the excited wave function is so small that the Rabi oscil- 
lations are hardly noticeable: within a Rabi period, most 
of the wave function has tunneled in the free region. The 
overlap plot of Fig. ||b confirms these statements, show- 
ing a strong decay over a few Rabi oscillations. Fig.||a 
illustrates how the switching of the current in the free 
region can be tuned by varying the width of the bar- 
rier. If we were to further reduce the width of the barrier 
(say, b = 0.5), the Rabi oscillations would cease to be no- 
ticeable: the width of the quasi level becomes large and 
the (on resonance) matrix element is strongly reduced 
in magnitude. Consequently, the activation mechanism 
provided by the dipole potential ceases to be efficient. 

We note that it is also possible to regulate the magni- 
tude of the leakage current by varying the off resonance 
mismatch Sui . According to Eq. (Q), changing the mis- 
match also changes the Rabi frequency. This is illus- 
trated in Fig. |J for the same coupling strength, we com- 
pare the integrated density for mismatchs Sui/qj = 0.01% 
(practically on resonance), Suj/uj = 5%, and Slu/uj = 
10%. Away from resonance, the ground state couples 
with continuum states neighboring the quasi level posi- 
tion, which have smaller matrix elements. Consequently, 
the leakage current is smaller than on resonance. More- 
over, as Su: is increased, the first term in the square root 
of Eq. ([I]) becomes important, and the period of the Rabi 
oscillations is decreased. 

We conclude with a brief discussion of the feasibility of 
this device. Our main concern is whether phonon relax- 
ation processes could induce transitions from the quasi 
level to the ground state before the excited states can 
tunnel out of the barrier. We proceed to convert the pa- 
rameters of the problem into dimension full units. If we 



2 



choose the height of the barrier to be of the order of leV, 
and the effective mass m* = 0.067m e for GaAs, the width 
of the well corresponds to a length a ~ 4.8 x 10~ 7 cm. The 
excitation frequency is u> ~ 0.8 x 10 15 s _1 , and the con- 
stant electric field associated with the dipole potential 
is E ~ 3 x 10 5 V/cro. For GaAs, typical optical phonon 
relaxation rates range from r p h ~ 10 _10 s to 10~ 7 which 
is very large compared to all other characteristic length 
scales of the problem. We therefore do not expect any 
significant changes in the qualitative behavior of the de- 
vice from these processes. 

In summary, we have proposed a switching device 
based on Rabi oscillations. The device is composed of 
a well region with a bound state and a quasi-level. Upon 
activation by a microwave field, a particle in the ground 
state can be activated to the quasi level, and subse- 
quently tunnel the free region adjacent to the well. The 
leakage current can be controlled in several ways: by 
varying the intensity of the driving field, by adjusting 
the width or height of the tunneling barrier, or by tun- 
ing the resonance mismatch to a desired value. We also 
note that this device could in principle also be used "in 
reverse" to measure the energy of a particle incident on 
the barrier and well region, from the free region. The in- 
cident particle tunnels in the well, and in the event that 
its energy corresponds to the quasi-level energy, it will 
get absorbed in the ground state. If the energy of the in- 
cident does not meet this resonance condition, the wave 
function should tunnel back in the free region without 
noticeable absorption in the ground state. 



FIG. 2. a) Integrated density in the "free" region as a 
function of time, for several values of the coupling constant: 
e = 0.1 (full line), e = 0.05 (dashed line), e = 0.2 (dotted 
line). The time is indicated by the number of infinitesimal 
time steps St = 0.0125 x (2n/u>) required for the evolution, 
b) Overlap with the ground state as a function of time, with 
the same input parameters. 

FIG. 3. a) Integrated density at fixed coupling strength, 
but for different barrier widths: b = 2 (full line), b — 1.5 
(dashed line), 6 = 1. (dotted line), b = 0.75 (dashed dotted 
line), b) overlap with the ground state, same parameters. 

FIG. 4. Integrated density, as a function of time for res- 
onance mismatch parameters 8lo/uj — 0.01% (solid line), 
Slo/uj = 5% (dotted line), and 5ui /u> = 10% (dashed line) 
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FIG. 1. Potential landscape: a well region is separated by a 
potential barrier from the "free region". The well depth and 
barrier height are adjusted so that there is only one bound 
state in the well, and there is only one quasi-level below the 
barrier 
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